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THE DETECTION OF SILICATE EMISSION FROM QUASARS AT 10 AND 18 MICRONS 
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ABSTRACT 

We report the spectroscopic detection of silicate emission at 10 and 18 /im in five PG quasars, the first detec- 
tion of these two features in galaxies outside the Local Group. This finding is consistent with the unification 
model for Active Galactic Nuclei (AGNs), which predicts that an AGN torus seen pole-on should show a sil- 
icate emission feature in the mid-infrared. The strengths of the detected silicate emission features range from 
0.12 to 1.25 times the continuum at 10 fim and from 0.20 to 0.79 times the continuum at 18 /Ltm. The silicate 
grain temperatures inferred from the ratio of 18-to-10/im silicate features under the assumption of optically 
thin emission range from 140 to 220 K. 

Subject headings: galaxies: active — quasars: emission lines — galaxies: ISM — infrared: galaxies 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are broadly classified in two 
types. Type 1 AGNs display broad hydrogen emission lines in 
the optical, whi le Type 2 A GNs do not. The AGN unification 
model (e.g. Antonuc cil 19931) ties these two types together. An 
optically and geometrically thick dusty torus surrounds a cen- 
tral black hole, accretion disk and broad-line emission region. 
Sources viewed face-on are recognized as Type 1, while edge- 
on sources are Type 2. Observations at many w avelengths 
support this scenario (e.g. Antonucci & Miller 1985), but until 
now there was a lack of constraining data in the mid-infrared. 

Silicate dust is a major component of the interstellar 
medium in the Milky Way and other galaxies. It produces 
two spectral features in the infrared, a "10 fim feature" which 
arises from a Si-O stretching mode, an d an "18 /nm feature" 
from an O-Si-O bending mode (e.g. iKnacke & Thomson! 
119731) . Emission or absorption from silicate dust is a domi- 
nant feature of the infrared spectrum produced by most mass- 
losing evolved stars, interstellar clouds, Hll regions, and cir- 
cumstellar disks around young stellar objects (YSOs). Sil- 
icate absorptio n is also d etected in dusty starburst galax- 
ies ( Rigopoulou et^y[199 l 9j) i _ ultraluminous infr ared galaxies 
(UI JR C.s: iC.enzel et al JI1 998t iTran et al]l200"Tl iSnoon et.alJ 
120041 [ Armus et al. 2004), and a sample of mid-infrared de- 
tected, optically invisible, high-luminosity galaxies with red- 
shifts of 1.7 < z < 2.8 (Houcket all2005l) . 

Type 2 AGNs show silicate dust in absorption as expected, 
but the mid-infrared spectroscopi c data for Type 1 AGNs 
have been scarce. Recently, Siebenmorgen et al. (2005) have 
presented spectra of two quasars with just the 10 /im feature 
in emission. While radiative transfer models clearly p redict 
silicates to be in emission for Type 1 sources (Laor & Draine 
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JT9931 iPier & Krolikl 119921 Il993l iGranato & Danesel Il994l 
Efstathiou & Rowan-Robinson ] Il995t INenkova et al] 
20001 1200 It 120021 Ivan Bemmel & Dullemondl 120031 
Dullemond & van Bemmel 2005), the absence of the 
silicate emission features at 10 an d 1 8 um has led to mo dels 
witho ut strong silicate emi s sion. ILaor & Draine ( 199 J) and 
van Bemmel & Dullemond (2003) have suggested that larger 
grains dominate the grain-size distribution. INenkova et alJ 
(2002), instead, proposed that if the torus is clumpy, the 
10 ^m s ilicate emission can be s ufficie ntly suppressed, 
although Dullemond & van Bemmel (2005) have challenged 
this suggestion. 

Here, we report the first spectroscopic detection of both 
silicate features at high signal to noise in five PG quasars. 
Quasars are the high luminosity counterparts of Seyfert 1 
galaxies, showing the same broad emission lines that define 
Type 1 AGNs. 

2. OBSERVATIONS AND DATA REDUCTION 

The five sources presented in this letter (PG 0804+761, 
PG 1211+143, PG 1351+640, I Zw 1 [=PG 0050+124] and 
3C 273 [=PG 1226+023]) were selected from a sample of 
12 AGNs based upon their prominant silicate emission fea- 
tures at both 10 and 18 fim. The observations were part of the 
Guaranteed Time Ob servation (GTO) pr ogram of the Infrared 
Spectrogra ph dRS~) 9 l|Houck et al.l2 004) on the Spitzer Space 
Telescope ( Werner et al. 2004). Table^lists the basic prop- 
erties of these targets, along with their observation dates and 
on-source in tegration times. All of them have broad H/3 lines 
iBoroson & Greenl ll992) and therefore are classified as Type 
1 AGNs. Two of them, 3C 273 an d PG 1211 are radio-loud 
quasars and the rest are radio-quiet dSanders et alJl!989i) . All 
five quasars are very luminous, as c an be seen from the bolo- 
metric luminosities in Table[0(from Sander s~t al.ll989l) . 

The observations were made with the Short-Low (SL) and 
Long-Low (LL) modules of the IRS. The spectra were ex- 
tracted from the flatfielded images provided by the Spitzer 
Science Center (pipeline version SI 1.0.2). The images were 
background-subtracted by differencing the two SL apertures 
and for LL, by differencing the two nod positions. Spectra 

9 The IRS was a collaborative venture between Cornell University and Ball 
Aerospace Corporation funded by NASA through the Jet Propulsion Labora- 
tory and the Ames Research Center 
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TABLE 1 
Properties of Sources 





PG 1351 


IZw 1 


3C 273 


PG 1211 


PG 0804 


AOR key 


3760640 


3761920 


4978176 


3760896 


9074944 


Date observed 


2004/4/17 


2004/1/7 


2004/1/6 


2004/1/7 


2004/3/1 


Integration time (s) 


408 


224 


336 


408 


408 


Redshift 


0.0881 


0.0611 


0.158 


0.0808 


0.0999 


Lum. distance (Mpc) a 


398 


271 


749 


363 


455 


log(L bo i/L ) 


12.14 


12.20 


13.44 


12.44 


12.56 


FlSum (Jy) b 


0.21 


0.62 


0.36 


0.21 


0.13 



"assuming H =71 kms~' Mpc" 1 , fi M =0.27, Oa=0.73, n K =0 
b rest frame 
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FIG . 1 . — The 5-35 /im IRS low resolution spectra of five PG quasars show- 
ing silicates in emission. The pixel-to-pixel uncertainty in the data is less than 
3%. Fringing effects with amplitudes of up to 5% are visible in LL order 1 
(see text). For each spectrum, three choices for the silicate-graphite contin- 
uum are indicated by the dashed line and grey shaded area. The spectra have 
been scaled for plotting purposes, with scaling factors of: 8.0 for PG 1351, 
1.0 for IZwl, 0.77 for 3C 273 and PG 121 1 and 0.50 for PG0804. 



were then extracted and calibrated using the IRS standard star 
HR 6348 for SL and the stars HR 6348, HD 166780, and HD 
17351 1 for LL dSloan et all2005l) . 

After extraction the orders were stitched, requiring order- 
to-order scaling adjustments of less than 5%. Finally, the 
stitched spectra were scaled to match the observed IRAS 
fluxes at 12 and 25 /im or IRS blue peak-up flux at 16 /mi. 
For 3C 273, we did not attempt any scaling, since the source 
varies by up to 0.4 magnitude s at 10.6 pm on time scales of 
up to several years (Neugebau er et alll9 99). The scaling fac- 
tors we applied average 6%, with the largest scaling factor of 
16%. The observed sources are point sources for both IRAS 
and IRS, therefore scaling to IRAS flux is justified. Figure 
[^presents the calibrated spectra. Table^lists the 15 /mi rest 
frame fluxes. 

In some spectra, most notably in IZwl, residual fringing, 
which are due to minor pointing errors, appears between 21 
and 30 /mi (the 1st order of the LL module). We have not 
defringed the spectra. Residual artifacts also appear in some 
spectra in SL order 1 (7.5-14 /mi). These take the form of a 
roughly sinusoidal deviation from the actual spectrum and are 
most obvious in the spectrum of 3C 273 at 8-10/mi, where 
they have a maximum strength of ^3% of the total flux. Be- 
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FIG. 2. — Continuum-subtracted silicate emission spectra of the five PG 
quasars. The spectra have been normalized and offset from each other for 
plotting purposes. Fringing effects are visible in LL order 1 (see text). The 
shaded areas indicate the uncertainty in the silicate emission spectra as intro- 
duced by the uncertainty in the observed spectrum and choice of continuum. 
Vertical dotted lines indicate the rest wavelengths of emission features dis- 
cussed in the text. 

cause 3C 273 has the weakest silicate emission feature in our 
sample, this artifact is even more pronounced in figures show- 
ing the continuum-subtracted spectra. We estimate uncertain- 
ties by comparing the spectra in the two nod positions, and 
these artifacts generally manifest themselves as a disagree- 
ment between the nod positions. The uncertainties only be- 
come apparant in Figures [2] and |3]shortward of 11 /mi (dark 
shaded area). None of these artifacts change the conclusions 
reported below. 

3. ANALYSIS 

The spectra in Figure are dominated by two broad emis- 
sion features, which we identify as the Si-O stretching mode 
and the O-Si-0 bending mode in silicate dust grains. The 
strength of the features is best appreciated by comparing the 
peak fluxes to the underlying silicate-graphite "continuum". 
For PG 1351, I Zw 1 and PG 1211, we define the contin- 
uum using a spline interpolation between the 5-8 spec- 
trum and a point at the extreme red end of the IRS spectral 
range, where a change in slope indicates the end of the 18 /im 
silicate feature. For 3C 273 and PG 0804, the spectral struc- 
ture at the red end of the IRS spectral range does not show 
a similar change in slope. This may indicate that dust emis- 
sion is still contributing at these wavelengths. Therefore, we 
have assumed that the continuum at the red end runs some- 
what below the observed spectrum. Figure^plots the result- 
ing spline-interpolated continua (dashed lines), while Figure 
|2]shows the resulting continuum-subtracted spectra. 

Comparison of the silicate emission spectra in Figure|2]re- 
veals that the ratio of the 18 and 10 /im features varies signif- 
icantly within the sample. It is highest for 3C 273 and lowest 
for PG 0804 (see Table |2}. Also the onset and center of the 
10 /mi features shift significantly from one source to the next, 
3C 273 and PG 1211 seem to appear at longer wavelengths. 
Interestingly, these two sources happen to be the only radio 
loud quasars in our sample. In order to test whether these 
differences can be attributed to our choice of continuum, for 
each source we compared the silicate emission spectra for two 
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FIG. 3. — The opacity spectra for the five PG quasars compared to the 
scaled synthetic silicate opacity curve of Weinaartner & Draine 1 200 1 ) and 
ILi & Dr aine 1 2001). The opacity spectra from the quasars are derived us- 
ing the adopted silicate grain temperature and scaled and offset for plotting 
purposes. The shaded areas indicate the uncertainty in the silicate profile as 
introduced by the uncertainty in the observed spectrum and choice of contin- 
uum. 



TABLE 2 
Measured quantities 



PG1351 I Zwl 3C273 PG1211 PG 0804 



10/im-feature-to-cont. ratio 11 


1.25 


0.38 


0.12 


0.55 


0.60 


18/tm-feature-to-cont. ratio 11 


0.79 


0.20 


0.36 


0.56 


0.56 


18/im/ 10/tm feature ratio 11 


1.75 


1.31 


3.34 


1.89 


1.08 


F(30/tm)/F(5.5/Ltm) 


11.0 


6.86 


2.53 


2.99 


1.95 


Lsil/Lbol 


0.087 


0.040 


0.009 


0.023 


0.021 


10/18/tm cont. color temp (K) b 


250 


270 


320 


310 


340 


Silicate temperature (K) 


200 


220 


140 


180 


210 



"evaluated at the peak of the feature 

b the color temperature is computed from the adopted continuum at 10 and 
18 /im 

other choices of underlying silicate-graphite continuum: one 
above and one below the adopted continuum, enclosing the 
shaded areas in Figure ^ As illustrated by the shaded areas 
in Figure[2] this results in an uncertainty of up to 30% in the 
peak flux of the 18 /im feature and in the ratio of the 18 to 
10 /im features. In contrast, the onset and center of the 10 /im 
features are not affected by different choices of continuum. 

The presence of the silicate features at 10 and 18 /im allows 
us to estimate the mean tempeature of the emitting grains by 
performing an unweighted \ 2 fit to the continuum-subtracted 
spectrum. This temperature is only physically meaningful un- 
der the assumption that the feature emitting region is optically 
thin. Using the continuum-subtracte d silicate opacity (see 
top pa nel i n Figure 13 derived from [Weingartner & Draine 
(2001) and ILi & Draind ( 120011) . we force the product of a 
blackbody and the silicate feature opacities to fit the observed 
continuum-subtracted silicate features. This reduces the prob- 
lem to a single free parameter: the grain temperature, and 
avoids concerns related to a possible nonthermal contribution 



to the continuum. The resulting silicate grain temperatures 
(T si] ) range between 140 K for 3C 273 and 220 K for I Zw 
1, see Table [2] These temperatures will change by as much 
as 15 K depending on the choice for the underlying silicate- 
graphite continuum. Uncertainties in the fitting process may 
amount to as much as another 10K. 

Using the derived silicate grain temperature (T S ji), it is pos- 
sible to compute the silicate opacity spectrum by dividing the 
observed silicate emission spectrum by the blackbody spec- 
trum with temperature T s ji. Figure [3] presents the resulting 
opacity spectra for our sample. 

The top panel in Figure|3]shows the synthetic silicate opac- 
ity spectrum after subtraction of the silicate-graphite contin- 
uum component. The lower panels show the opacity spec- 
tra for the five quasars overlaid on the synthetic profile. 
The quasar profiles generally show good agreement at wave- 
lengths beyond 1 /mi. Below 1 /im, only PG 1351 and I Zw 
1 match the synthetic profile well. In the other three spectra, 
the blue wing of the 10 /im appears weak, especially in 3C 
273 and PG 1211. Apart from the artifacts in the 8-10 /im 
region of the spectrum of 3C 273 (indicated by the large error 
bars in Figure|2j see §2) and the presence of atomic lines and 
PAH emission in some spectra, the 10 /im band has a smooth 
appearance, showing no sign of a departure from amorphous 
grain structure. Evidence for some crystalline grains do ap- 
pear at longer wavelengths, most notably in the spectrum of 
PG 1351 at 23 /im. 

The weakness of the blue wing of the 10 /im silicate opac- 
ity profile in three of our sources appears to be significant and 
cannot be explained by the uncertainties in the observed spec- 
trum (see the error bars in Figure[3}. Grain size and compo- 
sition, geometry, optically thick radiative transfer or a combi- 
nation of these factors may explain the observed weakness. If 
our optically thin assumption is invalid, then the temperatures 
we have inferred cannot be interpreted as real dust temper- 
atures. Jaffe et al. (2004) reported a similar deficiency in the 
blue wing of the lOum absorption feature in NGC1068, which 
they attribute to a different silicate grain composition. 

A close inspection of the spectra in Figures [2 and [2] re- 
veals numerous other features. Those visible in at least 
one spectrum in the sample include the PAH emission 
features at 6.2, 7.7 and 11.2 /mi, and emission lines at 
10.5 /tm ([Siv]), 12.8 /im ([Nell]), 14.3 /zm ([Nev]), 15.6 /tm 
([Nem]), 18.7 /im ([Sill]), and 25.9 /mi ([Oiv]). 

4. COMPARISON WITH MODELS 

Direct comparisons of our observations to the predictions 
of published AGN radiative transfer models are difficult given 
the limited resolution of model spectra in the 5^-0 /tm range 
of the IRS. Additionally, the model output is not tuned to 
be compared to the observed properties presented in Ta- 
ble|21 However, we can measure the featur e-to-continuum 
ratio a t 10 /mi from the figures pre s ented by lLaor & Draind 
(119931): iGranato & Danesd (11994 : iNenkova et aTn 2002): 
Dullemond & van Bemmel ( 2005). For the optically thin dust 
models of lLaor & Draind (119931) the values range from 0. 1 
for models where the grain-size distribution is dominated by 
larger grains to 1 .5 for models with a stan dard grain-size dis- 
tribution. In the face-on cases modeled by Granat o~& Danesel 
(1994), the values range from near zero for very thick and 
compact configurations to 1.0 for more exten ded torus mod- 
els. T hes e results are similar to those by INenkova et al.l 
( 120021) and lDullemond & van Bemmell (120051) . but their high- 
est 10 /tm feature-to-continuum ratios are around 0.7 and 0.4, 
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respectively. Given the observed range of 0.12 to 1.25 within 
our limited sample (Table |3, the model predictions for the 
10 /mi feature-to-continuum ratio are in good agreement. A 
similar comparison for the 18 /mi feature-to-continuum ratio 
could not be undertaken, since the features are not as apparent 
as the 10 /mi features in the AGN models. 

5. DISCUSSION AND CONCLUSION 

We have detected the 10 and 18 /im silicate emission fea- 
tures in five PG quasars — the first detection of both silicate 
features in emission in galaxies outside the Local Group and 
in AGNs in particular. Given the existence of published mid- 
infrared photometric and spectroscopic da ta for these quasars , 
our finding may appear surprising at first. Roche et al. ( 1991 ) 
obtained 8-13 /im spectra of 3C 273 and I Zw 1, but the low 
contrast of the silicate feature to the strong continnum and the 
lack of coverage outside the N-band window prevented them 
from identifyin g silicate e mission. The PHT-S spectrometer 
on ISO (Lemke et al. 1996) obtained spectra of 3C 273, 1 Zw 
1 and PG 0804 (Rigopoulou et al. 1999), but the limited wave- 
length coverage of PHT-S (only to 1 1.6 /mi) and the redshifts 
of the quasars allowed only the blue side of the 10 /mi feature 
to be visible. Because of these difficulties, silicate emission 
had remained unidentified in these five sources until now. 

Our detection of silicate emission in Type 1 AGNs pro vides 
support for the unified AGN model (Antonucci 1993J. An 
AGN torus viewed pole-on should offer an unobstructed line 
of sight to the torus dust, the surface of which is heated to 
temperatures of several hundred to a thousand Kelvin by the 
radiation from the central engine. Emission from this hot sur- 
face should produce a silicate feature in emission under most 
circumstances. 

Assuming a simple optically thin geometry, we find that the 



silicate grain temperature in our sources ranges from 140 to 
220 K, which is well below the dust sublimation temperature. 
Detailed comparison of the derived opacity profiles for our 
quasars indicates that our simple model cannot satisfactorally 
explain the absence of a blue wing in the 10 /im feature for 
three of our quasars. This likely points to a more complex 
geometry, involving radiative transfer within a (partially) op- 
tically thick environment - probably an AGN torus. It may 
also point to differences in the grain size or composition. Our 
derived silicate grain temperatures should, therefore, only be 
used as a characterization of the observed spectra, not as a 
means to infer the actual temperature and geometry of the 
AGN. 

From our limited sample of PG quasars (three in addition to 
the sample discussed in this Letter) and a few Seyfert 1 galax- 
ies, it is clear that 10 /im silicate emission is not a ubiquitous 
feature in mid-infared spectra of Type 1 AGNs. Only 5 out 
of 12 continuum-dominated AGNs show this feature promi- 
nently. Interestingly, the quasars which show a 10 /im sili- 
cate emission feature all have bolometric luminosities higher 
than 10 12 Lq, while the other Type 1 AGNs are less lumi- 
nous. Models which can explain the presence of the 10 /im 
silicate emission feature in some type 1 AGNs will also need 
to explain their absence in others. 
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